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Abstract 

We discuss the early evolution of ultrarelativistic heavy-ion collisions within a multi- 
fluid dynamical model. In particular, we show that due to the finite mean-free path 
of the particles compression shock waves are smeared out considerably as compared to 
the one-fluid limit. Also, the maximal energy density of the baryons is much lower. 
We discuss the time scale of kinetic equilibration of the baryons in the central region 
and its relevance for directed flow. Finally, thermal emission of direct photons from 
the fluid of produced particles is calculated within the three-fluid model and two other 
simple expansion models. It is shown that the transverse momentum and rapidity 
spectra of photons give clue to the cooling law and the early rapidity distribution of 
the photon source. 
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1 Introduction 



We discuss a three fluid hydrodynamical model for ultrarelativistic heavy-ion collisions in the 
energy range from BNL-AGS to CERN-SPS. The three fluids are introduced to separate the 
various rapidity regions observed in high-energy pp-collisions [1] . There it was found that the 
baryon charge essentially remains close to projectile resp. target rapidity and the energy loss 
due to particle production is transfered to midrapidity. Assuming that the initial (!) stage 
of ultrarelativistic heavy-ion reactions is essentially an incoherent superposition of binary 
iVW-collisions leads to a different picture of the early reaction dynamics as compared to 
one-fluid hydrodynamics. The projectile and target nucleons (as well as the newly produced 
particles) do not thermalize instantaneously (in the first interaction) but have to undergo 
several scatterings until kinetic equilibrium is eventually established. 

If the above-mentioned picture of the compressional stage of the heavy-ion reaction is 
correct, it is natural to introduce three fluids corresponding to the nucleons of the projectile 
and target (fluids one and two), and to the particles produced around midrapidity (fluid 3). 
The fluids have to be coupled via local friction forces leading to energy- and momentum 
exchange. The third fluid is, of course, absent initially and is produced in the course of the 
collision due to binary collisions between the nucleons of projectile and target. Presently, 
the rescattering of the produced particles with the nucleons is not taken into account. The 
interactions between fluids one and two are derived from data on iViV-collisions. We employ 
the parametrization of ref. [2] . A detailed discussion of our model as well as results on baryon 
stopping, kinetic equilibration, directed and radial baryon flow at AGS were presented in 
refs. [3]. In contrast to previous multi-fluid models [4], in our model 

1. three fluids, corresponding to projectile, target, and produced particles 

2. local interactions between the fluids which are based on iViV-interactions 

3. (3 + 1) -dimensional propagation of all three fluids in space-time 
are implemented. 

2 Compression shocks 

In (ideal) one-fluid hydrodynamics the projectile and target nucleons stop completely (in 
the center-of-momentum frame) when they touch. Therefore, two shock discontinuities are 
created [5] that propagate outwards if the velocity of the nucleons is supersonic and if matter 
is thermodynamically normal, cf. ref. [6] and references therein. Matter between the shock 
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Figure 1: Energy density profile of the projectile fluid along the beam axis (in the center 
of the reaction, i.e. x = y = 0) in one-fluid (left) and three-fluid (right) hydrodynamics as a 
function of CMS-time. 

waves is completely at rest (in the CMS) and reaches the maximal possible energy and baryon 
density that is compatible with local energy-momentum and baryon number conservation 
[7, 6]. In fig. 1 we show that the energy density of the shocked projectile matter and the 
velocity of the shock wave in the three-dimensional calculation agree reasonably well with 
the values obtained in the semi-analytical one- dimensional shock model (indicated by the 
box) [7, 6] with our equation-of-state, which treats the baryonic fluids as a non-relativistic 
ideal gas with compression energy [5]: 

2 

p(e,n) = -(e- E c n) +p c . (1) 
For the compression energy, we employ the ansatz 

E c = 7^— (n-n ) 2 + m N + W ,n ^0.1Qf m - 3 , (2) 
so that the compressional pressure p c is 

dE c k c , 2 2\ (o\ 
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However, one also observes that the full-step SHASTA [8] overshoots the correct value for 
the energy density of the compressed matter at the shockfront. This could probably be cured 
using a half-step treatment [6], which however requires considerably higher computing times 
and thus can not be easily implemented in multi-dimensional, multi-fluid models. Once the 
shock wave has reached the back-side of the projectile nucleus, a rarefaction wave starts to 
propagate inwards with the velocity of sound, leading to expansion and cooling [9]. This 
rarefaction wave can be clearly seen in the lower left panel. 

In the three-fluid model, the energy density profile looks qualitatively different, cf. right 
column of fig. 1. One observes that the projectile nucleons cross the symmetry plane z — 0, 
i.e. their mean-free path in the target matter is non-vanishing (in one-fluid hydrodynamics 
it follows from symmetry arguments that the baryon and energy flow through the plane 
z = vanish). Consequently, no real discontinuity appears, the energy density profile is 
smooth (within the numerical accuracy). The shock wave (if there is any) is smeared out 
considerably, which is due to the finite mean-free path (which in turn is linked to the total 
NN cross section). Also, the smeared-out shock wave obviously leads to less compression 
and heating of the shocked matter as compared to the discontinuity in the one-fluid limit. 
The entropy production in the two scenarios remains to be investigated. 

3 Kinetic equilibration of projectile and target nucle- 
ons 

One important issue in ultrarelativistic heavy-ion collisions is whether local thermal equi- 
librium is established via collisions between the particles. To study this, we have calcu- 
lated the ratio of the sum of the individual pressures of the projectile and target fluids, 
Pi(t, x) + p 2 (t,x), to the corresponding equilibrium pressure, p eq (t,x) (and similarly for the 
baryon densities). This ratio was then averaged over the volume where projectile and target 
overlap. This ratios approach unity if the relative velocity (at the given space-time point) 
between the projectile and target fluids becomes comparable to the thermal (resp. Fermi-) 
velocities of the particles within each fluid (the exact definition is given in ref. [3]). Thus, 
in one-fluid hydrodynamics it is by definition always equal to one. In contrast, in the three- 
fluid model (cf. fig. 2) ((pi +P2)/Peq) is close to zero in the beginning and reaches unity 
only after t^ M ~ IRp^j^cM = 1-5 fm/c. From this point on local kinetic equilibrium is 
established and the one-fluid limit is valid for the subsequent expansion. Although this time 
may look very short, it is significant for observables which are sensitive to the very early 
reaction dynamics, e.g. directed flow and "hard" thermal photons (or dileptons). Due to the 
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3f-hydro: Pb(160GeV)+Pb, b=3fm 
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Figure 2: Ratio of the sum of projectile and target pressure (and baryon density) to the 
equilibrium value as a function of CMS-time; Pb + Pfr-collisions (6 = fm) at CERN-SPS. 

kinetic non-equilibrium in the early stage, the conversion of initial longitudinal momentum 
into transverse momentum is weaker than in the one-fluid limit. 

4 Directed Flow at SPS 

We already discussed in ref. [3] that in this model kinetic equilibrium between projectile 
and target nucleons in the central region occurs on the same time scale as is relevant for 
the onset of directed nucleon flow. As a consequence, in the three-fluid model the directed 
nucleon flow is considerably lower and in particular less sensitive to the equation-of-state of 
the nucleon fluids than in the one-fluid limit. In fig. 3 we show our result for (pf r /N)(y) 
for Pb + Pfr-collisions at CERN-SPS (in calculating this quantity we have neglected thermal 
smearing of the nucleon momenta; for a definition see ref. [3]). The maximum is on the 
order of 40 — 50 MeV and considerably lower than at AGS (due to the fact that the time 
scale is shorter by roughly a factor of three). Quantitative comparisons to experimental data 
(directed flow was recently discovered in such reactions [10]) may reveal whether our present 
model produces enough directed flow or if the contribution of the pressure of the produced 
particles (as mentioned above, rescattering of produced particles with the projectile and 
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Figure 3: Directed nucleoli flow {pf r /N)(y) in Pfr + Pfr-collisions (6 = 3 fm) at CERN-SPS. 



target nucleons is presently neglected in our model) is essential. The produced particles are 
closer in rapidity to the projectile nucleons than the target nucleons (and vice versa), and 
thus might thermalize faster and exhibit significant pressure on them. Directed nucleon flow 
at SPS could thus be more sensitive to the equation-of-state of the produced particles than 
to that of the nucleons themselves. This will be investigated in future work. 



5 Thermal photon emission from the third fluid 

In our model, in Pb + Pfr-collisions at CERN-SPS the third fluid locally reaches energy 
densities up to pa 10 GeV/fm 3 (for comparison: in the one-fluid limit and for an equation- 
of-state appropriate for an ultrarelativistic ideal gas, p — e/3, the central region reaches an 
energy density of e/e = ^Icm ~ 3, where e Pd 0.15 GeV / fm 3 denotes the energy density of 
nuclear matter in the ground state). We therefore employ an equation-of-state for an ideal 
QGP (described within the MIT bagmodel) above T c = 160 MeV. Below T c , we assume 
that n, rj, p and u are the most abundant particles. The two equations-of-state are matched 
by Gibbs conditions for phase equilibrium, thus leading to a first-order phase transition. 

The temperature of the third fluid in the central cell is depicted in fig. 4. We also compare 
to two other hydrodynamical models for the expansion of the midrapidity fluid. In scaling 
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Figure 4: Temperature in the central region as a function of time in various hydrodynamical 
models, Pb + Pfe-collisions (b = fm) at CERN-SPS. 



hydrodynamics one assumes that the longitudinal flow velocity is v z = z/t [11] (independent 
of time), whereas the initial condition for the (one-dimensional) Landau expansion is v — 0. 
In these latter models we assumed an initial temperature of 300 MeV (for the boostinvariant 
expansion an initial time of t = 0.22 fm was employed). The maximum temperature 
in the three-fluid model is determined by the coupling terms between fluids one and two 
(and, of course, by the equations-of-state). The produced particles cool fastest in scaling 
hydrodynamics, while in the case of a Landau expansion the temperature is constant until 
the rarefaction waves reach the center. In this latter case one finds the longest-lived mixed 
phase, due to the fact that we assumed purely one-dimensional (longitudinal) expansion. 

The various cooling laws reflect in the transverse momentum distribution of direct pho- 
tons [12], which can be produced in the QGP- (by annihilation of quarks and antiquarks 
and Compton-scattering) and hadronic phase (mainly by 7T — p scattering) [13]. The faster 
the cooling of the photon source, the steeper the slope of the photon spectrum. Fitting 
the spectrum in the region 2 GeV < k T < 3 GeV by an exponential distribution we find 
"T" = 260 MeV in the three-fluid model and for the Landau expansion, and "T" = 210 MeV 
for the scaling expansion. 

The rapidity distribution of thermal photons with transverse momenta much larger than 
the maximal temperature of the photon source (e.g. kr = 2 GeV at SPS) directly reflects 
the rapidity distribution of the photon source at early times [14], cf. fig. 6. This is due to the 
fact that "hard" photons do not thermalize, in contrast to e.g. pions. Therefore, in the three- 
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Figure 5: Transverse momentum distribution of direct photons at midrapidity (i/cm = 0) 
within various hydrodynamical models, Pb + Pfo-collisions (6 = fm) at CERN-SPS. 

fluid model and in the Landau-expansion case the photon rapidity distribution is strongly 
peaked around midrapidity and is not proportional to the (squared) rapidity distribution of 
the pions at freeze-out, as in scaling hydrodynamics. 

6 Conclusions 

We discussed the compression of projectile and target and showed that the shock disconti- 
nuities occuring in the one-fluid limit are smeared out considerably in the three-fluid model. 
The energy density of the baryons is less than in the one-fluid limit. Most of the energy 
loss due to stopping goes into production of new particles, leading to the creation of a very 
hot (temperatures up to 300 MeV at CERN-SPS) third fluid. We compared the cooling of 
this third fluid in the three-fluid model and two other simple expansion models, longitudinal 
Landau expansion and longitudinally boostinvariant (plus cylindrically symmetric radial) 
expansion, and showed how the various cooling laws reflect in the transverse momentum 
distribution of thermal photons. The rapidity distribution of thermal photons with trans- 
verse momenta much larger than the maximal temperature measures the rapidity spread of 
the photon source at early times (before acceleration leads to broadening). Finally, in the 
three-fluid model the time scale for kinetic equilibration of the nucleons in the central region 
is found to be on the order of the Lorentz-contracted nuclear diameter, 2R/^ CM . This leads 
to a less efficient conversion of initial longitudinal momentum into transverse momentum as 
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Figure 6: Rapidity distribution of direct photons with kx = 2 GeV calculated within various 
hydrodynamical models, Pb + Pfr-collisions (6 = fm) at CERN-SPS. 

compared to the one-fluid limit. As a consequence, the directed flow, which is also created 
in the early stage of the reaction, is less than in the one-fluid limit. 
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